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Abstract. Wc consider the Landau-Lifshitz equations of ferromagnetism (in- 
cluding the harmonic map heat-flow and Schrodingcr flow as special cases) for 
degree m equivariant maps from to §P. If m > 3, we prove that near- minimal 
energy solutions converge to a harmonic map as i ^ oo {asymptotic stability), 
extending previous work [11] down to degree to = 3. Due to slow spatial decay 
of the harmonic map components, a new approach is needed for to = 3, involving 
(among other tools) a "normal form" for the parameter dynamics, and the 2D 
radial double-endpoint Strichartz estimate for Schrodinger operators with suffi- 
ciently repulsive potentials (which may be of some independent interest). When 
m = 2 this asymptotic stability may fail: in the case of heat-flow with a fur- 
ther symmetry restriction, wc show that more exotic asymptotics arc possible, 
including infinite-time concentration (blow-up), and even "eternal oscillation". 
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1. Introduction and Results 

The Landau- Lifshitz (sometimes Landau- Lifshitz- Gilbert) equation describing the 
dynamics of an 2D isotropic ferro magnet is (eg. [12]) 

= ai(An + |Vnpn) + X An, ai > 0, 02 G M (1.1) 

where the magnetization vector u = u(t,x) = {ui,U2,U2) is a 3-vector with normal- 
ized length, so can be considered a map into the 2-sphere 

u : [0, T) X ^ §2 ._ £ ^3 I I -| ^ ^^ 2) 

The special case a2 = of (1.1) is the very well-studied harmonic map heat-flow into 
S^, while the special case ai = is known as the Schro dinger flow (or Schrddinger 
map) equation, the geometric generalization of the linear Schrodinger equation for 
maps into the Kahler manifold 

In order to exhibit the simple geometry of (1.1) more clearly, we introduce, for 
u G S^, the tangent space 

Ts§2 := = e I M ■ f = 0} (1.3) 

to the sphere at u. For any vector G M^, we define two operations on vectors: 

r:=vx, P'^:=-rr. (1.4) 

For u G S^, P" projects vectors orthogonally onto T^jS^, while J" is a 7r/2 rotation 
(complex structure) on T^§^. Denoting 

a = ai + ia2 e C, (1.5) 

the Landau-Lifshitz equation (1.1) may be written 

ut = P:Au, := aiP" + 02 J" (1.6) 

The energy associated to (1.1) is simply the Dirichlet functional 

S{u) = 1-1 \Vu\^dx (1.7) 

and (1.6) formally yields the energy identity 

S{u{t)) + 2ai [ [ \P'^Au{s,x)\^dxds = S{u{0)) (1.8) 
Jo 

implying, in particular, energy non-increase if ai > 0, and energy conservation if 
ai = (Schrodinger map). 

To a finite-energy map -u : ^ is associated the degree 

deg(n) := 7- / u^^ ■ J^'Ux^dx. (1.9) 

If lim|2;|^oo {{(x) exists (which will be the case below), we may identify u with a 
map §^ — *■ S^, and if the map is smooth, deg(u) is the usual Brouwer degree (in 
particular, an integer). It follows immediately from expression (1.9) that the energy 
is bounded from below by the degree: 

£{u) = 11 \u,, - J"m,,|2 + 47rdeg(M) > 47rdeg(M), (1.10) 
2 Jr-2 
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and equality here is achieved exactly at harmonic maps solving the first-order equa- 
tions 



ru., (1.11) 



which, in stereographic coordinates 



§23^ „ !fl±l!fa e C u {00} (1.12) 

are the Cauchy-Riemann equations, and the solutions are rational functions. These 
harmonic maps are critical points of the energy S and, in particular, static solutions 
of the Landau-Lifshitz equation (1.1). 

In this paper we specialize to the class of m-equivariant maps, for some m G Z"^: 

u{t,x) = e"''^^v{t,r), {7: [0,T) X [0,cx)) ^ §2 (1.13) 

with notations 

i?:=/ = A;x, A? =(0,0,1), (1.14) 

and polar coordinates 

xi + 1x2 = re*^. (1.15) 
In terms of the radial profile map v = {vi,V2,V3), the energy is 

£(n) = 7r^" (^\v^.\^ + !I^(vl + vl)^rdr. (1.16) 

Finite energy implies v is continuous in r and limr^o v = ±k, lim^^oo v = ±k 
(see [10] for details). We force non-trivial topology by working in the class of maps 

:= {u = e'^^^t;(r) | S{u) < 00, v{0) = -k, v{oo) = k}. (1.17) 

It is easy to check that the degree of such maps is m: 

deg = ^- (1.18) 

The harmonic maps saturating inequality (1.10) which also lie in E.^ are those 
corresponding to (3z"^ {(5 G = C\{0}) in stereographic coordinates (1.12). In the 
representation §^ C M^, the harmonic map corresponding to is given by 

e"^^MO, (/ii,0,/i3), /ii = , /i3 = . (1-19) 

The full two-dimensional family of m-equivariant harmonic maps in is then 
generated by rotation and scaling, so for s > and a G M, we denote 

^i = m\ogs + ia, h[fi] h' = h{r/s). (1.20) 

The harmonic map e"^^^h[fi\ corresponds under stereographic projection to e~^z"^. 

We are concerned here with basic global properties of solutions of the Landau- 
Lifshitz equations (1.1), especially the possible formation of singularities, and the 
long-time asymptotics. 

For finite-energy solutions of (1.1) in 2 space dimensions, finite-time singularity 
formation is only known to occur in the case of the 1-equivariant harmonic map 
heat-flow (a2 = 0) - the first such result [4] was for the problem on a disk with 



4 



STEPHEN GUSTAFSON, KENJI NAKANISHI, AND TAI-PENG TSAI 



Dirichlet boundary conditions (this was extended to Si on in [9]). Examples of 
finite-time blow-up for different target manifolds (not the physical case S^) are also 
known (eg. [20]). 

For the Schrodinger case (oi = 0), [5] showed that small-energy solutions remain 
regular. In the present setting, the energy is not small - indeed by (1.10) and (1.18), 

£ k„ > 47rm. (1.21) 

A self-similar blow-up solution, which however carries infinite energy, is constructed 
in [6]. 

In the recent works [11, 9, 8], it was shown that when m > 4, solutions of (1.1) 
in Em with near minimal energy {£{u) ~ Airm) are globally regular, and converge 
asymptotically to a member e"^^^h[fj] of the harmonic map family. In particular, 
the harmonic maps are asymptotically stable. The analysis there fails to extend to 
m < 3, due to the slower spatial decay of -^h[fi] (a point which we hope to clarify 
below). With a new approach, we can now handle the case m = 3 as well: 

Theorem 1.1. Let m > 3, a = ai + ia2 G C\{0}, and ai > 0. Then there 
exists 6 > such that for any u{0,x) G with S{u{0)) < imn + 5^, we have a 
unique global solution u E C([0, oo);Sm) of (1.1), satisfyingVil E Lj ^^^([0, cxd); L^). 
Moreover, for some fi E C we have 

\\u{t) - e™^^/I[/i]||Lgo + axE{u{t) - e'"^^K[/i]) ^0 as t ^ oo. (1.22) 

In short, every solution with energy close to the minimum converges to one of 
the harmonic maps uniformly in x as t ^ oo. Even for the higher degrees m > 4, 
this result is stronger than the previous ones [11, 9, 8], where the convergence was 
given only in time average.^ Note that in the dissipative case [a\ > 0), solutions 
converge to a harmonic map also in the energy norm, while this is impossible for 
the conservative Schrodinger flow (ai = 0). 

The analysis for the case m = 2 seems trickier still, and we have results only in 
special case of the harmonic map heat-fiow (02 = 0) with the further restriction that 
the image of the solution remain on a great circle: ^2 = (a condition preserved 
by the evolution only for the heat-fiow). These results show, in particular, that 
the strong asymptotic stability result of Theorem 1.1 for m > 3 is no longer valid; 
instead, more exotic asymptotics are possible, including infinite-time concentration 
(blow-up) and "eternal oscillation" : 

Theorem 1.2. Let m = 2 and a > 0. Then there exists 6 > such that for any 
u{0,x) = e'^^^v{0,r) G S2 with £{u{f))) < 811 + 6^, and V2{0,r) = 0, we have a unique 
global solution u G C([0,oo);S2) satisfying Vu G ^^^^^([0, 00); L^). Moreover, for 
some continuously differentiable s : [0, 00) (0, 00) we have 

\\u{t) - e™^^7l(r/s(t))||igo +£(M(t) - e™^^/I(r/s(t))) ^0 as t 00. (1.23) 

In addition, we have the following asymptotic formula for s{t): 

(1 + 0(1)) log(s(t)) ^ 2 /--^ M^^^ ^ (1.24) 

TT Ji r 



^The statements in the previous papers do not follow directly from Theorem 1.1, but are implied 
by the proof in this paper. 
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where as t ^ oo, o(l) — and Oc(l) converges to some finite value. In particular 
there are initial data yielding each of the following types of asymptotic behavior: 

(1) s{t) 3 e (0,oo). 

(2) s{t) ^ 0. 

(3) sit) oo. 

(4) = liminf < limsups(t) < oo. 

(5) < liminf s{t) < limsup s{t) = oo. 

(6) = liminf < limsups(t) = oo. 

Estimate (1.23) shows that these solutions do converge asymptotically to the 
family of harmonic maps. However, the evolution along this family, described by 
the parameter s{t), does not necessarily approach a particular map in S2 (although 
it might - case (1)). The solution may in fact converge pointwise (but not uniformly) 
to a constant map ±k (which has zero energy, zero degree, and lies outside S2) as 
in (2)- (3) (this is infinite-time blow-up or concentration), or it may asymptotically 
"oscillate" along the harmonic map family, as in (4)- (6). 

Note that the above classification (l)-(6) is stable against initial "local" pertur- 
bation. Namely, if two initial data v^{0) and w^(0) satisfy 



the corresponding solutions have the same asymptotic type among (l)-(6). More 
precisely, the difference of their scaling parameters converges in (0, 00). The point 
is that the energy just barely fails to control the above integral. 

In particular, the oscillatory behavior in (4)-(6) is driven solely by the distribution 
around spatial infinity. In fact, if we replace the domain by the disk D = {x & 
I \x\ < 1} with the same symmetry restriction with m = 2 and the same boundary 
conditions v{t, 0) = —k and v{t, 1) = k, then it is known [1, 7] that all the solutions 
behave like (2), namely they concentrate at x = as t ^ 00. Also, if we replace the 
domain by S^, then we can rather easily show in the dissipative case ai > that 
the solution converges to one harmonic map for all m G M, by the argument in this 
paper, or even those in the previous papers. We state the result on in Appendix 
A with a sketch of the proof. 

As for other target spaces, we should mention an example in [19, Section 5], 
which is "eternal winding" by heat flow around a compact 1-parameter family of 
harmonic maps from to x with an artificial warped metric. Our result has 
the following comparative advantages: 

(1) The geometric setting is very simple and physically natural. 

(2) The asymptotic formula is explicit in terms of the initial data. 

In addition, our analysis works in the same way in the dissipative (ai > 0) and the 
dispersive (ai = 0) cases. We need 02 = in Theorem 1.2 only because the angular 
parameter a{t) gets beyond our control (hence we remove it by the constraint), but 
the rest of our arguments could work in the general case.^ 

^We will use the parameter convergence in the proof of Theorem 1.1 in the dispersive case ai = 
to fix our linearized operator. However it is possible to treat the linearized operator even with 
non-convergent parameter and oi = 0, if we assume one more regularity on the initial data. We 
do not pursue it here since the wild behavior of a{t) prevents us from using it. 




(1.25) 
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1.1. The main difficulty and the main idea. The standard approach for asymp- 
totic stabihty is to decompose the solution into a leading part with finite dimensional 
parameters varying in time, and the rest decaying in time either by dissipation or 
by dispersion. In our context, we want to decompose the solution in the form 

v{t) = h[ii{t)]+v{t) (1.26) 

such that the remainder v{t) decays, and the parameter yu(t) G C converges as t ^ oo 
(at least for Theorem 1.1). In favorable cases (the higher m, in our context), we can 
choose ii{t) such that all secular modes for v{t) are absorbed into the time evolution 
of the main part /i[/x(t)]. This means that the kernel of the linearized operator for 
v{t) is spanned by the parameter derivatives of and hence we can put that 
component of dtv{t) into dth[^{t)\. This is good both for v{t) and yu(t), because 

(1) v{t) will be free from secular modes, and so we can expect it to decay by 
dissipation or dispersion, at least at the linearized level. 

(2) The decomposition is preserved by the linearized equation. Hence /i(t) is 
affected by v only superlinearly, i.e. at most in quadratic terms. 

In particular, if we can get decay of v in time, then /i(t) becomes integrable in 
time, and so converges as t ^ oo. This is indeed the case for m > 3. 

However, the above naive argument does not take into account the space-time be- 
havior of each component. The problem comes from the fact that the decomposition 
and the decay estimate must be implemented in different function spaces, and they 
may be incompatible if the eigenf unctions decay too slowly at the spatial infinity. 

In fact, the parameter derivative of is given by 

dh[^i] = Ke'^^lihl, 0, -hl)dfii + (0, 1, 0)df^i2\ (1-27) 

and hence the eigenfunctions are 0(r~™) for r oo, i.e. slower for lower m. On 
the other hand, the spatial decay property in the function space for the time decay 
estimate is essentially determined by the invariance of our problem under the scaling 

v{t,x) ^ v{X'^t,Xx), (1.28) 

which maps solutions into solutions, preserving the energy. If we want decay in 
time (so that we can integrate quadratic terms in /i), then a function space with the 
right scaling is given by 

v/r G L^L^. (1.29) 

To preserve such norms in x under the orthogonal projection, the eigenfunction 
must be in the dual space, for which m > 3 is necessary. Indeed, this is the essential 
reason for the restriction m > 4 in the previous works [10, 11, 9]. We emphasize 
that the above difficulty is common for the dissipative and dispersive cases, since 
they share the same scaling property. That is, the dissipation does not help with 
this issue, even though it gives us more flexibility in the form of decay estimates. 
The main novelty of the present approach is the non-orthogonal decomposition 

Ll = (K) © {^n^, (1.30) 

where ip^{r) is smooth and supported away from r = and from r = oo, so that 
the (non-orthogonal) projection may preserve the decay estimates. This is good 
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for the remainder v, but not for the parameter /i — the decomposition is no longer 
preserved by the hnearized evolution, since they have no particular relation. This 
implies that we get a new error term in fi{t) which is linear in v(t) (see Section 6). 
This contribution is handled by including it in a sort of "normal form" for the 
dynamics of the parameters explained in Section 7. In particular, it is this 
new term which drives the non-trivial dynamics for the m = 2 heat-flow given in 
Theorem 1.2. 

For the purely dispersive (Schrodinger map) case, one tool we use should be of 
some independent interest: the 2D radial "double-endpoint Strichartz estimate" for 
Schrodinger operators with sufficiently "repulsive" potentials (in the absence of a 
potential, the estimate is false). The proof is given in Section 10.2. 

1.2. Organization of the paper. In Section 2, we use the "generalized Hasimoto 
transform" to derive the main equation used to obtain time-decay estimates of the 
remainder term. Section 3 gives the details of the solution decomposition described 
above, and addresses the inversion of the Hasimoto transform. The estimates for 
going back and forth between the different coordinate systems (the "Hasimoto" one 
of Section 2 and the decomposition of Section 3) are given in Section 4. Section 5 
is devoted to establishing the time-decay (dispersive if Oi = 0, diffusive if Oi > 0) 
of the remainder term, using energy-, Strichartz-, and scattering-type estimates. 
The dynamics of the parameters n{t) are derived and estimated in Section 6. The 
leading term in the equation for /i is not integrable in time, and so Section 7 gives 
an integration by parts in time to identify (and estimate) a kind of "normal form" 
correction to /i(t), whose time derivative is integrable. At this stage, the proof of 
Theorem 1.1 for m > 3 is complete. A more subtle estimate of an error term for 
m = 3 is done in Section 8, completing the proof in that case. Finally, in Section 9, 
the normal form correction is analyzed in the case m = 2, a2 = 0, V2 = 0, in order to 
prove Theorem 1.2. Proofs of certain linear estimates (including the double-endpoint 
Strichartz) are relegated to Section 10. Appendix A states the analogous theorems 
for domain and sketches the proofs. 

1.3. Some further notation. We distinguish inner products in and C by 

3 

a-6 = ^^afc6fc, a ° 6 = Re a Re 6 -|- Ima Im6. (1-31) 

k=l 

Both will be used for C'^ vectors too. The inner-product is denoted 

{f\9)= [ f{x)W)dx. (1-32) 

while (/,(?) just denotes a pair of functions. For any radial function /(r) and any 
parameter s > 0, we denote rescaled functions by 

r(r) := /(r/.), /^(r) := f{r/s)s-\ (1.33) 

We denote the Fourier transform on by JF, and, for radial functions, the Fourier- 
Bessel transform of order m by JF^: 

i^fm = ^[ f{x)e-^^<dx, {Tjf{p)= r Jmirp)f{r)rdr, (1.34) 
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where Jm is the Bessel function of order m. For m G Z we have 

■Jm{r) = ^ r e'^'-'''''^'de, ^[/(r)e™^] = z"^(^„/)e^"^^ (1.35) 

We denote the Laplacian A^. on the subspace spanned by d-dimensional spherical 
harmonics of order m by 

A^™^ := dl + {d- l)r~^dr - m{m + d- 2)r~^ (1.36) 

Finally, the space is the dyadic version of U\rdr) defined by the norm 

ll/IU. = ||||/(r){2^- < r < 2^-+i}|U.(,,,)|L,(^). (1.37) 

2. Generalized Hasimoto transform 

First we recall from the previous papers [10, 9, 8] the equation for the remainder 
part, which is written in terms of a derivative vanishing exactly on the harmonic 
maps, and so independent of the decomposition. The equation was originally derived 
in [5] in the case of small energy solutions (hence with no harmonic map component), 
and called there the generalized Hasimoto transform. 

Under the m-equivariance assumption (1.13), the Landau-Lifshitz equation (1.6) 
is equivalent to the following reduced equation for v{r,t): 



(2.1) 



Define the operator d^j on vector- valued functions by 

777 ^ 

■= dr rR. (2.2) 

r 

Since for any vector 6, J^'Rh = k{v -b) — {k ■ v)b, we have J^R = —v^ on the tangent 
space T^S^ = v-^. For future use, we denote the corresponding operator on scalar 
functions by 

■.= dr + ^. (2.3) 

r 

Then equation (2.1) can be factored as 

V, = -P:Dld,v, (2.4) 

where 

D ■= P^dP^ (2.5) 

will always denote a covariant derivative (which acts on T^S^-valued functions), and 
* denotes the adjoint in L^(M^). Denote the right-most factor in (2.4) by 

m 

w:=d^^v = Vr P^'k. (2.6) 

r 

Then (2.4) becomes Vt = —P^Dlw, and applying to both sides yields 

Dtw = -P^D^Dtw. (2.7) 



STABILITY AND OSCILLATION OF HARMONIC MAPS 9 

Now we rewrite the equation for w by choosing an appropriate orthonormal frame 
field on T^S^, reahzed in C'^. Let e = e(t,r) satisfy 

Reeev^, I Reel = 1, Ime = J'^Ree. (2.8) 

Let S, T be real scalar, and let q, v be complex scalar, defined by 

w = goe, P'^k = uoe, Dte = -iSe, Dre = -iTe. (2.9) 

Then we have the general curvature relation 

[Dr, Dt]e = i{Tt - Sr)e = -i det {v Vr Vt) e. (2.10) 

Using the equation (2.4) for v, we get 

Sr-Tt = {w + -P^k) ■ PlDlw. (2.11) 

Now we fix e by imposing 

D^e = 0, e(r = oo) = (1, i, 0). (2.12) 

(The unique existence of such e will be guaranteed by Lemma 4.1.) Then (2.11) 
yields 

Sr = {q + —v)°{iaL*;q). (2.13) 

A key observation is that in the Schrodinger (non-dissipative) case a = i, we can 
pull out the derivative on q: Sr = {dr + ^){^\q\'^ + ° q), and so 

S = Q-r2Q* (3:=i|# + !!l..,= iH' + I!i^^l = (2,14) 

Jr r 2 r 2 r 

The evolution equation (2.7) for w yields our equation for q: 

poo 

I 771 

{dt - iS)q = -aL^Llq, S = - (g + -i^) ° {iaLlq)dr. (2.15) 

J r ^ 

This is the basic equation used to establish diffusive (ai > 0) or dispersive (ai = 0) 
decay estimates. The operator acting on q can be expanded as 

L^Ll = d;dr + - ^ + , ^ + — W3. (2.16) 

3. Decomposition and orthogonality 

Next we introduce coordinates for the decomposition of the original map 

v = h[^i]+v, (3.1) 

or more precisely for the remainder v, and a localized orthogonality condition which 
determines the decomposition. The choice of coordinates is the same as in the 
previous works [11, 9, 8], while the decomposition itself is different. 

For each harmonic map profile fi = m log s+ia, we introduce an orthonormal 
frame field 

f = f[^]:=e"«(-/I^x/ + ^/). (3.2) 
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on the tangent space T^j^jS^, such that the parameter derivative of h[fi] is given by 

dh[iJ,] = hldi^oi. (3.3) 
We express the difference from the harmonic map in this frame by 

z:=v-i. (3.4) 

In other words P'^^f^^v = z°i,OTV = z°i + "yh[fi], where we denote 



7:= vr^N'-i = -o(kr). (3.5) 

As explained in the introduction, the orthogonahty condition in the previous works 

{z I hi) = (3.6) 

would not work for m < 3 due to the slow decay of hf for r ^ 00. Hence instead 
we determine the parameter /i by imposing localized orthogonality 

{z I if') = 0, ip' = ip{r/s), (3.7) 

with some smooth localized function (f{r) G C^((0, 00); M), satisfying {hi \ (f) = 1. 
The fact that e^^^/if/i] solves (1.11) means that 

d^h = 0, (3.8) 

and so we have 



Hence 



TTl TtX 

w = d^v = Vr ^ {h'^^v + -1)3 -u) = Uv H v-^v. (3.9) 



777 ?77 

Uz = Uv-i + v-ir = w-i- —V3Z + —hlj. (3.10) 



In order to estimate z hj w (or equivalently g), we introduce a right inverse of the 
operator = dr + ^h^, defined by 

poo pr 

R^g.-h\{r) / / h\{r"Y^g{r")dr"i^^{r')h{{rydr' (3.11) 



^0 Jr' 

Then we have 

L'R'^g = g, R'^L' g = g - hl{g \ , (3.12) 
hence = {L"^)"^ on {f'^)'^- Moreover we have the following uniform bounds 
Lemma 3.1. For all p G [l,C)o] and \6\ < m, we have 

WRtpdWrOL^ ^ IIV^IIr-eLiJIfi'llre+iLi, (3.13) 

where the norm is defined in (1.37). Moreover, the condition on (f is optimal in 
the following sense: if f > 0, then ip G r^^L^, is necessary for R^ to be bounded 
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We give a proof in Section 10. Note that the above bounds are scahng invariant: 
denoting Dsf := f{r/s), we have 



i?; = sD,rId:\ = s~'DsL'd;\ 

We can combine the estimates of the Lemma with the embedding 

r'^m C r'^m ^ l-e, = --e2, pi> P2, qi < q2- 



The above lemma is used as follows. First note that the orthogonality | z 
implies that z = R^L^z because of (3.12). For the energy norm, we choose 6 = 
and p = 2 in Lemma 3.1. Then 

Since — dr\ < 1/r, we further obtain 



i?; : r'L' ^ r'^X, {\e\ < m), 
where the space X is defined by the norm 

\\4x ■■= \\z/r\\L2 + WZrWLl- 

The Sobolev embedding X C L°° is trivial by Schwarz: 

Iklligo < ||2;/r||i2||z^||i2. 

Hence we get by using (3.10), 

Ikllx < ||^''2;||i2 < \\q\\L2 + \\z\\l^\\z\\x- 
For Lf estimates of z, we use Lemma 3.1 with 6=1 and p = oo. Then we 

\\z/r\\L^ < \\ip\\r-iLl\\L'z\\r2Li < \\L'z\\l^, 

for any p G [1, oo], and so by using (3.10), 

Ik/^ II L?L- < Ik II L?L- + IklU- |k||L?L^^- 



3.14) 
3.15) 

) = o 



3.16) 
3.17) 



3.18) 

3.19) 

3.20) 
have 
3.21) 

3.22) 



If we were to use hi instead of ip, then we would need m > 3 for the Strichartz- 
type bound (3.22), and m > 2 for the energy bound (3.16), by the last statement of 
the lemma. 

4. Coordinate change 

Before beginning the estimates for the evolution, we establish in this section the 
bi-Lipschitz correspondence between the different coordinate systems: v and (fJ^jq), 
including unique existence of the decomposition. It is valid for any map in our class 
Em with energy close to the ground states. 

For that purpose, we need to translate between the different frames e and f . At 
each point (t, r), we define = f ® e G GL^{C), a real-linear map C ^ C, by 

Mz:=i-{e^z). (4.1) 

Its transpose tM = e (g) f , defined by ^Aiz = e ■ (f ° 2), is the adjoint in the sense 
that (Aiz) °w = z° (^Aiw). For any b, c, d G C'^ we have 

b-(cod) = (Reb-c)od + (Imb-c)od. (4.2) 
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— Mv (g) e. 
r 



(4.3) 
(4.4) 



Since f{oo) = e~^°'e{oo), and f ± we have 

M{oo) = e~*", Mr = i (S) Br + ir e = -i v{e ■ u 
Then e can be recovered from Ai by 

e = pM/^le + (h[fj] . e)h[fi] = 'Mi - (1 + -fy\Mz)h[fi] 
provided that I7I < 1. We further introduce some spaces with (pseudo-)norms. 
S^{v) = £(e'"^^w(r)), = min(| Re/i|, 1) + dist{lmfi, 27rZ), 

\\4x = \\z/r\\L2 + \\Zr\\Ll, W^Wy = WMrhHdr) + \\M\\l 

J^miS) = Mr) : [0,00] 
L\6) = {q{r):[0,oc)- 
The metric on C is defined such that 

\\h[^i'] - h[^']\\x - IIM/] - h[fi']\\L^ - 1/ - (4.6) 



52 I v{0) = -k, v{oo) = k, Sm{v) < 4m7r + 5^}, 
: I \\q\\L2 < V26}, C = C/27rzZ. 



(4.5) 



Lemma 4.1. Let m G N and (f G Co(0, 00) satisfy {(p \ hi] 
(5 > such that the system of equations 



1. Then there exists 



v = h[ii\ + v, z = v-i\ii], (z I v^"") = 0, 7 = vl 

1)^6 = 0, e(oo) = (1, i, 0), 



r 



(4.7) 



defines a bijection from v G Sm(5) to (/i, g) G C x -^^(5), which is unique under 
the condition \\z\\l°° ^ S. v, z and e are also uniquely determined. Moreover, if 



e-') with j = 1,2 are such tuples given in this way, then we have 

,2 A^l A^2| 



\\v' - v^Wx + \\z' - z^Wx + \\e' - el^oo + \\M' - M'\\y 



(4.8) 



where M^ := P eK 

In particular, we have pointwise smallness 

II^IIl- ~ IklU- < 5 < 1, 

so that we can neglect higher order terms in z or 

Proof. We always assume (3.9), (3.4) and (3.5), which define the maps 

V ^ w = q° (-u, /i) H-i> ti ^ 2; 7, (i), fi) V, 

with the Lipschitz continuity 



(4.9) 



\x<\\z'-z'\\x-^\\v' 



v'Wx, 



\v'-v^x-\\v'-v'\\x\<\f^'-f^\. 



(4.10) 



(4.11) 



The energy can be written as 



2Srr,(v) 



m 



-Rv 



m 



-P^k 



Ll 



(4.12) 



\l, + 2m{vr I P'^k/r) = \\w\\l2 + A7i[vs{oo) - Vs{0)]. 
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Since ||P^'t/||x ^ £m{vy^'^, X C and \v\ = 1, the boundary conditions ^3(0) = —1 
and ^3(00) = 1 make sense in the energy norm. 

Next we consider a point orthogonahty. Let v G 2^(5). Since ^3(0) < < Vs{oo) 

and v^{r) is continuous, we have v{so) = e*"°^/i(l) for some /iq = mlogso + iao, so 
that V = /i[/io] + w is a decomposition satisfying {z \ ip^°) = if ip{r) = 6{r — 1). In 
this case v is recovered from {w, fio) by solving the ODE: 

L^z = wJ ■ f [H - -^3^2 + -/ii°7, ^(so) = 0, (4.13) 



or the equivalent integral equation 



z — 

^ ~ -^5(r-l) 



7T7 Tn 

w-i\ix^\--v^z^-h\'^ . (4.14) 



The uniform bound on R\(^^_^-^ can be localized onto any interval / 9 sq, because z 
is the solution of the above initial value problem. Hence we get, in the same way as 
in (3.16), 

IN||rL2nL-a) ~ II^IUlW + lklUg°(/)lkllrL2{7)- (4.15) 

Since z{sq) = and ||g||L2 < 5 ^ 1, we get by continuity in r for J — » (0, 00), 

IkllxnL- < IklUi < 5. (4.16) 

Thus every v G T,m{S) is close at least to some h[fio], and we have — ■if ^ X hj 
(4.11). Em(^) is a complete metric space with this distance. 

Now we take any ip G Cq{0, 00) satisfying [ip \ hi) = 1, and look for fi around /iq 
solving the orthogonality 

F(/i) := (tJ- f[/i] I ^^) = {v . f[/.] I ^^) = {z I ^^) = 0. (4.17) 

Its derivative in /i is given by 

dF = -{v-hlh[fj] I ip^)dfx - t{v- hli[fx] I ip^)da - iv-i[fA \ (rdr + 2)^^) — 

= -dfi - {v ■ h[fi] I hlif^)dfi - {i! ■ i[fi] I {rdr + 2)ip^ds/s + ih^ip^da) ^^"^^^ 
= —dfi + 0{6\dfi\). 
In particular we have 

\F{fio)\<S, |^(/io) = -/ + 0(5). (4.19) 

In addition, both F{fi) and df^F are Lipschitz in v. Therefore by the implicit map- 
ping theorem, if (5 > is small enough, there exists a unique G C for each v such 
that F{fi) = and |/i — /io| < 5, and v ^ fi is Lipschitz. Then 

\\4l^ ^ \W- MHIU? + l/io - /i| < 5 < 1, (4.20) 
and so by the same argument as for (4.16), we get ^ and in addition, 

Ik' - z'Wx < l/i' - /i'lc + 11^' - w^L^. (4.21) 
If we have two such /i = /ii, fi2 with Hz-' H^oo < 6, then 

l/ii - f^2\c ~ - MHIU- < 11^?- MHIU? + IK"^- ^NIU- < 5, (4.22) 
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and so the implicit mapping theorem imphes that fii = fi2- Thus we get a bijection 
V ^ (/i, w) with the Lipschitz continuity 

11^^' - v^Wx ~ l/i' - /i'lc + 11^' - w^hl- (4.23) 

For the frame field e, we consider the matrix = f e, together with the 
equivalent set of equations (4.3) and (4.4). Integrating (4.3) from r = oo, we get 

\\M - e-Wr < IKVrlUill^-'.llLi + \\v/r\\Li\\K/r\\,. < 6, 

\\M' - WWy < l/i' - /x'lc + 6\\v' - v^x + 6\\e' - e^^^ + \\v' - v^x, 

while (4.4) provides 

l|e' - eligo < \\M^ - M^Wl^ + |/ - fi% + \\z^ - z^\\l^. (4.25) 

Hence for fixed v G ^^(5) (and fi), we can get {Ai, e) G F x L°° by the contraction 
mapping principle for the system of (4.3) and (4.4). Moreover we get 

\\M' - M^Y + lle^ - elz.- < \\v' - v^x- (4.26) 

If (/i, g) G C X L'^{S) is given, we consider the system of equations (4.3), (4.4) and 

z = R 



Mq - —v^z + —Kl\ . (4.27) 

which is equivalent to the q equation in (4.7) under the orthogonality [z \ ip^ = 0. 
The last equation provides, through the uniform bound on R^, 



' <\f^'-f^'\c+\\q'-q'\\Li + S\\M'-M'\\Y + S\\z'-z'\\x. (4.28) 



\z - z Wx 



Combining this with (4.24) and (4.25), we get {z,M.,e) for any fixed (/x, g) by the 
contraction mapping, and moreover they satisfy 

\\z' - z^x + \\M' - M'Wy + \\e' - e'h- < |/ - fi% + \\q' - q^hi- (4.29) 

□ 

5. Decay estimates for the remainder 

In this section we derive dissipative or dispersive space-time estimates of the 
remainder v in terms of z, from the equation (2.15) for q. First by the smallness of 
z, we obtain from (3.20) and (3.22), 



X 



<Mt)\\L^<6, |k/r||L- < IklU-, (5.1) 
for all p G [1, oo]. Next we estimate the factor 5*, by using 

/oo 
fgdrhr <J2J2^'''\\f9\\LHr-2'^) ~ Wfah^ < WfhAlgiy- (5.2) 
jez k>j 

Then from the expression in (2.15) for S, we have 

\\smq < \\sit)i.^^Lr < ihhi + \\4rLi + mmui < ii^kiui. (5.3) 

In the dispersive case ai = 0, we avoid the derivative by using expression (2.14) 

WSmLl < \\Smr-^Lr < ikhl + MrLl + l)||g||L,- < khf' (^1 = 0). (5.4) 
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For the time decay estimates, we treat the dissipative and the dispersive cases sep- 
arately. 

5.1. Dissipative Lf estimate. Here we assume ai > 0. By the equation (2.15) of 
q, we have 

dtWqWh = -2ai\\Ltq\\l2, (5-5) 

hence 

||g|U-L2 + ||Ltg||^2i2 < ||g(0)||L2 ~ 5. (5.6) 
Since = I and -R^* : rL"^ by Lemma 3.1 and duahty, we have 

MrLl < ||^'*g||L2 < WLUhl + ||^^||L-||g||rL2- (5.7) 

Since the last term can be absorbed by (4.9) smallness of w, we get 

IklU < h/rhl + W*q\\Ll < \\L*MlI- (5.8) 
So by using the bound (3.17) on i?^, we obtain 

lkllrL?x < llg||L?x < WLkhi^ ^ lk(o)IUi ~ 5, (5.9) 

and also from (5.3), 

\\S\\^Ll<S. (5.10) 

5.2. Dissipative decay. Next we show the convergence q ^ as t —>■ oo, by 
comparing it with the free evolution. For T > 0, let 

g^:=g-e(*-^)'^^^'""^g(T). (5.11) 

Then we have 

qj - aAf = {^S - aV)q, q^T) = 0, (5.12) 
where the potential V{t,x) is given by 

2m(l — fq) m /r 1 ON 

V = ^ , + -W3. (5.13) 

Multiplying the equation with g^, we get the energy identity 

fa,{\\qJ\\l. + \\'^q^\\h)dt = Re [\-aVq + ^Sq \ q^)dt, (5.14) 
and hence by Schwarz, and using estimate (5.3) to put 5* G L'^L'^, 

^ k/rh^rl + ||g||L-L?|k||L?_X ^0 (T ^ oo), 



Hence ||g'(t)||L2 can not converge to a positive number, since e^* '^)°-^2 \{T) 
as t —>■ oo for all T > 0. Thus we obtain 

Mt)\\x<\\qmLi^O (t^oo). (5.16) 
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5.3. Dispersive Lf estimate. Next we consider the case Oi = (and 02 7^ 0). We 
set (with no loss of generahty) a = i. Since the energy identity provides only 
bound on q, we have to work with the Strichartz estimate in a perturbative way. 
Denoting H'^ := L'^L'^* , the equation of q is given by 

qt + iH'^'^k = Ni + N2, (5.17) 

where 

Ni := —2am— — - — g, N2:=iSq — 2am-^q — am—q, (5.18) 

and S is given by (2.14). We have 

m<\h,{s{t)/smMlr\ (5.19) 

and so 

< \\hMt)/<m\LTh\\LlLT- (5.20) 

Using (5.4), we have 

WSqWilLl < WSWilLihWi^.L^ < M\%L^- (5.21) 
The other terms in are bounded in L^L^ by 

Iklli^L- + lk/'^IL?L^lkllL?L- < h\\%L^- (5.22) 

Now we need the endpoint Strichartz estimate for with fixed scaling s\ 

Lemma 5.1. Let = L^L^* = —A^^^^^ + 2mr~^(l — h^) and m > 1. Then we 
have 



\„-iH''t,„\\ ^ 



II / e ^* ^^f{t')dt'\\L^LlnL^L^^\\f\\LjLl+L^Ll, 
J — 00 

uniformly for any fixed s > 0. 



(5.23) 



This Lemma will be proved in Section 10.2. Hence if | log(s(t)/s(0))| ^ 1 for all 
t, then we have 



hWhr Lin^L^ ^ lk(0)||L2 ~ 5, (5.24) 

and also from (5.4) 

\\S\\lIli<5- (5.25) 

5.4. Dispersive decay. Next we prove the following asymptotics of scattering type 
for q and z: 

e-ii^i"''"\{t) ^ 3q_^ in Ll z ^ in L~ {t ^ 00). (5.26) 
For the scattering of q, we further expand the equation 

qt-zAt-'\ = No + N2, (5.27) 
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where N2 is as in (5.18), and 

iVo := -2am^—-^q (5.28) 

Then the global Strichartz bound implies that 

\\No\\l^,l\{t,oc) ^ 0, \\N2\\lilI{t, 00) (5.29) 

as T — 00. By Strichartz (for once again, we get the scattering of q. 

For the vanishing of we use the inversion formula 

z = R'^g, g = Mq + r''^m{hl'y - vsz). (5.30) 

Since i?^ is bounded L°°, the latter two terms contribute at most with 

||-2||Lg° ||-2/r||L2 ^ ll^ll^^oo, hence we may drop them. Also we may replace q by 

. , , {m — 1) 

its asymptotic free solution := e 2 g^. Moreover we may approximate g+ 
by nicer functions. Hence we assume that q := e C^(0, 00). Then we may 

further replace the free solution with the stationary phase part: 



°°{t,r)=C^t-V/^''^ / J„_i(rp/(2t))e^''V(^*)g+(p)pdp 



(5.31) 

C^rV/^''^q{r/{2t)) + n, 



where the error is bounded by Plancherel 

||7^|U. ~ 11(1 - e^^'/(^*))g+|U. < t-irVlUi - 0. (5.32) 

Now that spatially local vanishing is clear (eg. it follows from \\R^Mq{t)\\rL°° ^ 
||q'°°(^)||l°° 0), we may extract the leading term of i?^ for large x. We assume 
that s{t) G and suppy?'* C (0, b) for a fixed b G (0, 00). Then for r > 6 we have 



iK9)ir) = 0(1) + f ^9{r")dr" 



(5.33) 



hl{r"y 

= 0(1) + / [r /r)^g{r')dr' as r — > 00. 
Jb 

Thus we are reduced to showing that 

Gx ■■= J\p/rrM{t,p)t-'e'P'/^''h{p/t)dp in (5.34) 

for any x ^ ^"^(0, 00). By partial integration on [p/t)e^P^ /^^^\ we have 
r^Gx = {i/2)[p^~^M{p)eP"'^''h{p/t)][ 

[{m-l)M{p)x{p/t)/p + M{p)x!{p/t)/t 

+ Mr{p)x{p/t)\p'^~'e^'"^^''^dp, (5.35) 

The right hand side is bounded by r™/t, using \x{.plt)\ ^ p/i for the first, second 
and fourth terms, |x'(p/^)l ^ 1 for the third, and M.r £ L'^L^{dr) for the fourth 
term. Thus we obtain ||2;(t)||j;,cx) — * 0. 
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6. Parameter evolution 

It remains to control the asymptotic behavior of the parameter fi{t) of the har- 
monic map part of the solution. Its evolution is determined by differentiating the 
localized orthogonality condition 

= dtiz I ip') = {vfi\ + {v-it\ + {z I dtif^, (6.1) 

and each term on the right is expanded by using 

Vt = Vt- h[^]t = -{aL*jq) o e - /i^/i f , 



it = -ih^ai - h{fih[fi], dtip" = --rdr'^". 

s 



(6.2) 



Plugging this into the above and then dividing it by s^, we get 

ll = -{MaLlq I - (hlfi^ I - (z I {^rdr - 1(12^)^^). (6.3) 

m 

The last two terms are bounded by 

lAllklU°°(llv5||Li + WrdMli^), (6.4) 

and so absorbed by the left hand side 
Since = \P^k\ ^ hi + \v\ and hence 

\vr\<\q\ + \z\/r + hl/r, (6.5) 

we get from (4.3), 

l-^rl < \qz\ + \z\yr + \zhl\/r. (6.6) 

The leading (first in the r.h.s) term in (6.3) can be estimated, using [JUa, L*.] = 
^Ara, as follows 

liMaL*^ I ^^)| < s-\\\MrU + ll-^IU^)lkll.Li 

(6.7) 

^ S (IklUi + \\z/r\\L2 + l)||g||rL2. 

Hence using that H^H^oo < \\z\\x < MIl^ < 5 < 1, we get 

II^^AIL? ^ h/^hl, ~ IklL^L^- (6.8) 

7. Partial integration for the parameter dynamics 

Now we want to integrate in t the right hand side of (6.3), which is not bounded 
in Lj. The key idea is to employ the q equation (2.15), by identifying a factor of 
L:ffLtq, through a partial integration in space. 

For the spatial integration, we first freeze the phase factor Ai. Since h[n\ = v = 
—k at r = 0, we have Ai{t, 0) = e''°, i.e. f (t, 0) = e*"e(t, 0) for some real a{t). Then 
Dti{t, 0) = ia'{t)f{t, 0) - iS{t, 0)f (t, 0), and so 

a'{t) = S{t,0) + a'{t). (7.1) 

We decompose 

M=e'^ + M, (7.2) 
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and rewrite the leading term of (6.3) as follows. Let c = ||/ii||^2- Since L^j = 
L* + mvs/r and L^hf = 0, we have 

{MaLtq I = ae'^iLtq \ + {MaL*jq \ 

= ae^° [{Ltq \ - ch,Y) + {mqv^/r \ c/i/)] (7.3) 
+ {Maq I L^Lp'^) + {Mraq \ ip^). 

The second term is bounded by Hgtr^^Hii < ||g/r||i2 ||2;/r^||i2, and the last two 
terms are bounded by 

\\q/r\\^{\\Mr/r\\L2 + \\M/r\\L^). (7.4) 
where the last factor is further bounded by using that Ai = aX r = 

\\M/r\\L^ < \\Mr/rh. < \\q/r\\Li + \\z/r'Ul < Ml^- (7.5) 
We further rewrite the remaining (main) term. By the definition of c, we have 

- ch, \ h) = 1 - c\\h,\\l, = 0, (7.6) 

and so we have 

ip^-chl = L^*R;*{ip'-chl), (7.7) 
where the operator was defined in (3.11). Let 

ij := R* {^ - chi) = ^r^— + 0(ri-3-) (r ^ oo), (7.8) 

V- m — 1 

where the asymptotic form easily follows from the fact that 

POO 

^{r) = -c{hi{r)rY^ I hi{r'fr'dr' (r > 1). (7.9) 



Then we have, by using equation (2.15) for g, 
{-aLtq\{^ - ch,Y) = {-aULtq \ ^ / s) 

= {qt - iSq I tp^/s) + {amq \ L^tJgr" VV'S), 
and, using (3.9), the last term is bounded by 

\m\ + \vlr\)r'%Ll < ih/rhl + Wz/r'^.f < \\q\\l^. (7.11) 
For m > 2, & L^, and so 

WiSq I r/s)\\Ll < WSWl^.lMl^l- < Shh^^L-, (7.12) 
either by (5.10) or (5.25). The last two terms of (6.3) are bounded in Ll by 

\\s(^\\L^\U/r\\^LT\\{r\v\ + r'\^r\Y\\Ll < MI^l^, (7.13) 

where we used (6.8) and (5.1). 
Thus we have obtained 

\\fi-e^Hqt\r/s)\\L}<S\\q\\L?L-. (7.14) 
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Integrating by parts in t, the leading term is rewritten as 

\qt I r/s) = dtie'^q \ ^1 s) - i{sa + sS{t, 0))(e*"g | ^P^) 
+ s{e^^q I {rdr + l)^l)^). 



(7.15) 



The last term can be bounded in L\ by using (6.8), 

\\sii\\L2\\{q I {rdr + < Ml^l^Uq \ (r9, + 1)^^)11^2. (7.16) 

If m = 2 or m > 3, then {rdr + l)"?/" ^ L^, and so the above is further bounded by 
II Q'll ^22,00- When m = 3, we need some extra effort to bound the last factor in - 
this is done in the next section. 

If m > 2, we have for the leading term 

|(e^"g|V^7.)|<||g|U.||V;|U2 <||g(0)|U2, (7.17) 

while for m = 2 this term can be infinite from the beginning. We will show in 
Section 9 that the time difference [(e*"g | ip'^/s)]^ can be controlled for finite t, but 
still may become unbounded as t — oo for some initial data. 

This also means that the second last term of (7.15) is beyond our control when 
m = 2, and so in this case we force it to vanish by making the assumptions 02 = 
and V2 = 0. For the other cases (m > 2), we should estimate S{t,0), for which we 
use in the Schrodinger case (a = i) that 

S = Q- -^dr, Q = -\q\' + -ws = 0(|gp + \z/r\^ + qhl/r), (7.18) 
Jr r 2 r 

since |u;3| < |g||i^| and |i^| ^ + \z\. Thus we get at each t, using (5.1), 

ll-^IU,- < \\q\\lf + s-'\\q\\L^. (7.19) 
Then the second term in (7.15) is bounded in 



g|li2..||(g I r/s)\\Lr + (II^«IIl? + lkllL?Lso)||(g | ^011 



(7.20) 



<M^L-{S+Uq\ij^)h.), 

where we used (6.8). If m > 3, then ip & and hence the last factor ||(g | V^^Hi.z 
is bounded by Ho'llia^^go ^ S. Its estimate for m = 3 is deferred to the next section. 
In the dissipative case ai > 0, we estimate simply by (2.15) at each t 

\\S\\l^ < ih/rhl + \\z/r'Ui + \\hl/r'h.)\\Llq\\L^., (7.21) 
and hence the second term in Lj is bounded by 

IklU^L^IkllL^x + (l|s/i||L2 + lkllL2x)ll(g I ^^)IIl2 

<ML^L^{S+\\{q\ij^)\\L.), 

where we used (6.8) and (5.9). 

Thus we have obtained all the necessary estimates to prove Theorem 1.1 when 
m > 3. Furthermore, its proof for m = 3 will be complete once we show 

Wiq \ + Wiq \ (rdr + l)iJ^)h^ < 6, (7.23) 

which will be done in Section 8. 



(7.22) 
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For Theorem 1.2, it remains to derive the asymptotic formula (1.24) from the 
leading term (g | ip^/s), and to show that all of the asymptotic behavior (l)-(6) can 
be realized by choice of the initial data u{0,x) - this is done in Section 9. 

8. Special estimates for m = 3 

In this section we finish the proof of Theorem 1.1 by showing (7.23). It suffices 
to estimate the leading term for r ^ oo: 

\\ir^\ir)\q)\\,.<5, (8.1) 

with X £ C°° satisfying x(r) = for r < 1 and x(r) = 1 for r > 2, since the rest 
decays at slowest 0(r~^) G L^, for which we can simply use q G L'^L'^. Once the 
above is proved, we can conclude that 

ll/^lUr ^ 1/^(0)1 + llglli-LinL^L-- (8.2) 

The boundedness of fi and the scattering of q imply that the "normal form" correc- 
tion (e*"g I ip^ / s) converges to zero, and so /i(t) is convergent as t — oo. 

To estimate (8.1), we use perturbation from the free evolution e"*^2 ; 

q - aA^q = N^ + iVa, (8-3) 

where A'^o and N2 are as in (5.28) and (5.18), satisfying 

iVo e r-2(r/.)-%2Lr, N, ^ L]LI. (8.4) 

For the contribution of A^2 as well as the initial data, we use the following estimate. 

Lemma 8.1. For any I > 0, anya G withRea > 0, and any functions g{r) , f (r) , 
and F{t,r) , we have 

U9\e'''^''f)huo,oo)<\\r'9\\LM\Ll, 

\\i9\ f e'^(*-^)^'V(.)rf.)|U.(^)<||r2^|U^.HF|Uj,.. ^^'^^ 

Proof. We start with the estimate for the free part. Let 'g = Tig and / = Tif . The 
above Li norm equals by Plancherel in space. 



e 



"OO 

"'*^7)|Il?(o,oo)~II / e-'^*'^G(a)da|U.(o,oo), (8.6) 







where we put 



G{a):=g{a')f{a^). (8.7) 
If ai > 0, then (8.6) is bounded by Minkowski 



/•oo nOQ 

< II / e-''^"'\G{a)\dah.^o,oo) < / e~''^-\\r'G{a/t)\\L.^,^^^da 

Jo Jo 

POD 

< 11^11^2(0,00) / a'^^^e'^^'^da < 1 1^1 1^2 (0,00) • 
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If Oi = 0, then 02 7^ and (8.6) is bounded by Plancherel in t, 

POO 

< II / e-*«^*'^G(a)rfs|L2(i,) ~ \\G\\l2^o,oo). (8.9) 
Jo 

Thus in both cases we obtain 

>(0 



e'^*^^7)IU?{o,oo) < l|G(a)|U.(M) < ||?||l^||/||l^ ~ ll?llL^II/llLi. (8.10) 

Then the first desired estimate follows from 

f°° f°° dr 

\g{p)\< \Ji{rp)\\g{r)\rdr<yg\\L^ \Ji{r)\ VgWi^, (8.11) 

Jo Jo 

since |Ji(r)| < min(r', r~^/^) for r > 0. 

By duality, the estimate on the Duhamel term is equivalent to 

II / Ks + t)e''''^^ g{x)ds\\L^Li<yg\\L^\\M\^. (8.12) 

which is equivalent to 

11/ A(t)e^*^^ <7(x)c?t|U. <||r2^|Ugo||A|L., (8.13) 

which is dual to the first estimate. □ 
For the potential part Nq, we transfer the equation to by n = r~^g and consider 

Ut-a/\fu = r-^No. (8.14) 
Then thanks to the decay of the potential, we have 

r-^iVo e ^2^^0/7(^6) ^ L^ijg-f (M^)^ (8.15) 

as long as s{t) is away from and 00. Then by the endpoint Strichartz or the energy 

estimate on M^, the corresponding Duhamel term is bounded in L^H^ ^''^(M^), and 

since |Vxr~^x('^)l ^ we have r~'^x ^ -^5/4 (J^^) -^3/2 (^^)- Thus to summarize, 
we have 

W-\ I q)\\^ < \\qm\Ll + IklLfL,- • (8.16) 

Completion of the proof of Theorem 1.1: Let initial data u{0) be specified 
as in Theorem 1.1. The existence of a unique local-in-time solution u{t) in the 
given spaces can be deduced by working in the (/i, q) variables (using the bijection 
of Lemma 4.1) and using estimates similar to those of Sections 5 and 6. The details 
are carried out in the Schrodinger case (a = i) in [11], and carry over to the general 
case in a straightforward way (in fact, there are well-established methods for energy- 
space local existence in the dissipative case, starting with the pioneering work [17] 
on the heat-flow). It follows from this local theory that the solution continues as 
long as fi{t) is bounded and q is bounded in Lf'Ll n L'^L^. 

For m > 3, the estimates of the previous four sections give the boundedness of q 
and n which ensure the solution is global, as well as the convergence of fi{t). The 
convergence to a harmonic map then follows from the estimates of Section 5. □ 
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9. Special estimates for m = 2, a > 0, v2 = 

Let m = 2 and (with no further loss of generahty) a = 1. By the bijective 
correspondence v <-> {fi,q), it is clear that W2 = is equivalent to /x, g G M. It 
remains to control the leading term for the parameter dynamics 

{q I r/s). (9.1) 

In particular, we will show that this can diverge to ±oo, or oscillate between them 
for certain initial data. 

First by the asymptotics for r ^ oo, we have ip + cr]~^ G L^, where we denote 

_i j r^^ (r > a) -i J ^"^ < ^) -i J ^ ""^ {a < r < b) 
'''^< " jo (r < a) |o (r > b) " |o (otherwise) ^ ' ' 

Hence we may replace s'^ip'^ by — cr~^ modulo o(l)L^. 

Next we want to replace q by the free solution g° := e*^2 'g(O). For that we use 
the following pointwise estimate to bound q — qo'- 

Lemma 9.1. Let Rea > and I > 1. Then for any function g{r) satisfying 
\g{r)\ < {r)~'^ , we have 

\e-t^^2^ g(r)\ <mm{l,l/r,r/t). (9.3) 
Proof. Let ai = Re a. By using the explicit kernel we have 

poo pTV 2 2 

e""'^' g{r) = - / e-" ^ '"^r^^^ +^^'g{q)qdedq, (9.4) 
and the integral in 6 can be rewritten by partial integration on e*^ as 



-Isinee"" -+'^^g{q)qdedq. (9.5) 

The double integral for |g — r|>r/2is bounded by using the second form by 

2 2 2 

/ ^e""^^^ min(g, l)dg < re~"^^ min(t^^/^t~^) < min(l,r/t, 1/r), (9.6) 
Jo ^ 

and that for |g — r|<r/2is bounded by using the first form by 

/ e-'''te-^''^demm{r,l)dq<t-h^/\t/r^Y/^mm{r,l) 
Jo J-n (9.7) 

< min(l, 1/r), 

and by the second form by < r^/t xl/r = r/t. □ 
The nonlinear part of g contributes as 



\q-q') = - (e^^'""''(*-*'V;(,V | V(t')g(tO)rft', (9-8) 

t/ 

where the potential term is given by 
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The contribution from the last two parts is estimated with the bound from the 
above Lemma, thus bounded by 

(IIW^Il?,, + U/rhOWq/rhl^ < \\q{0)\\%. (9.10) 

We need to be more careful to estimate the other term q{l — hD/r"^. First by 
Schwarz and the pointwise estimate, we have 

2 



J 



<\\q/r\\% l I mm{s{t)-\r-\T/{t-t')f{r/s{t'))-''^-dt' 

2m 



^0 



(9.11) 



r 



where we also used that |1 — h^{r)\ < (r) It suffices to bound the last double 
integral. Let t = t — t'. For < r < s(t)^, the r integral is bounded by 

/ --+ -J^-+ ^-<s{tr\l + \og{s{tf/r)\ (9.12) 



'T/s{t) ■^i'-; ' Js{t) 
hence its r integral is bounded by 



/ s(t)-2(l + log(s(t)Vr))dr = 1 + / |log0|d^<oo. 
Jo Jo 



For s(t)^ < r, the r integral is bounded by 



s{t') ^2 ^ p ^2^^^/)4m ^ ^^^/)2 



j-2 y, J s{t') r'^r^"^ r ^ t"^ 



(9.13) 



(9.14) 



and its r integral is bounded by square of 

\\s{t')/T\\Llis{t)\t) < ||s(t)/r||L2(,(j)2,t) + ||(s(t) - S(t - r))/r||i2(,(4)2^4) 

<i + rii=^iiL^^-'/'rf^<i. 

Thus we obtain 

\{r-4)< I ? - < lk/^llL?JII^3/rl,2^ + ||s|U2 + 1) < ||g(0)|U2, (9.16) 

namely we may replace q by the free solution g° in the leading asymptotic term. 
Furthermore, we can freeze the scaling parameter because 

<l[iog4r/Mi)<o(i)(l[iog4l + i)- 

Thus we obtain 

(l + o(l))[21og4 = -c[(r;(J)Jg°)]*+0(l), (t ^ oo) (9.18) 
where 0(1) is convergent. 
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The leading term is further rewritten in the Fourier space by using that 



oo 



•^i[-;(0)<](P) = P"W Mr)dr 

Js{o)p (9.19) 

= p-'Ms{0)p) = p<}/,(0) + s{0)R{s{0)p), 3R e Ll 
Let go '■= ^iQ'(O). By Plancherel we have 11^(0)11^2 = Hgollil and 

[-(^;(o)< I ?°)]o = ((1 - e-'-'Kl/sio) I 5)) + 0(1) 

= (-^<v.<iM0) 1^(0)) + 0(1) (9-20) 
= 27T q{0,r)dr + O{l). 

Js{0) 



2 _ ^-2^, 



Thus we obtain (using that c = ||/ii||j;^2 = tt 



l + o(l))[log4 = - [ q{0,r)dr + O{l), (9.21) 

71" Js(0) 



's(O) 

and the error term 0(1) converges to a finite value as t — *• oo. 

Completion of the proof of Theorem 1.2: As in the proof of Theorem 1.1, 
we now have all the estimates to conclude the solution is global (in particular, fi{t) 
remains finite by the above formula and estimates), and the convergence to the 
harmonic map family follows from the estimates of Section 5. It remains to consider 
the asymptotics of s(t). 

Since g(0) G does not require \q{0,r)\dr < oo, it is easy to make up 
q{0) G L^, for any given s(0) G (0, oo), such that the first term on the right of (9.21) 
attains arbitrarily given limsup > liminf G [—00,00] as t ^ 00. In particular, all 
of the asymptotic behaviors (l)-(6) in Theorem 1.2 can be realized by appropriate 
choices of (^'(0), s(0)), for which Lemma 4.1 ensures existence of corresponding initial 
data u(0) G S2. 

Using that f2 = 0, we can further rewrite the leading term in terms of v. Since 
e = (173, i, —Vi), we have 

q = we = ViVsr - ^s^ir + ~ = ^Pr + — , iy-22) 
where P is defined by -u = (cos P, 0, sin P). Hence we have 

(l+o(l))[log4 = - /''^M^dr + 0(1), (9.23) 

where 0(1) converges as t — 00. 
□ 



10. Proofs of the key linear estimates 
10.1. Uniform bound on the right inverse R^. 



(lo.i; 



:io.2) 
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Proof of Lemma 3.1. Let s = 1 and omit it. It suffices to prove 

\\R^g\\reL^ < \\(p\\j.-eLi\\g\\re+iLl^, 
ll-R^/llr-s-iL-^ ^ IIV^IIr-eL3„ll/llr-''L3„- 

From this we get by duality, 

ll-Rv'S'llrflLf ^ ||'y^||r-eLi„ ||5'||r«+iLi, 
ll-R^/llr-e-iLf ^ ||V^||r-eLi||/||r-9Li, 

and the bihnear complex interpolation covers the intermediate cases. 
It remains to prove (10.1). We rewrite the kernel of 

K9 = II ^^X(r, r', r")W)hi{ryg{r")dr"dT\ (10.3) 

where %('") is defined by 

{1 (r' < r" < r), 
-1 (r < r" < r'), (10.4) 
(otherwise). 

We decompose the double integral dyadically such that r ~ 2-' , r" ~ 2'^ and r' ~ 2\ 
and let 

Aj = 2'^^\R^g\\Lr^(^r^23), -Bj = 2^''+-^^-'||_R*/||ioo(r^2i), (10 5) 

= 2^'llV5||Li{r'~2')' 9k = 2^~^""^^''||fi'||Li(r^2'=)5 fk = 2^^'||/||Li(r~2fe)- 

For R^, we have 

J-1<A,<«+1 V-^w.w; 
l-l<k<j+l 

The sums over /c are bounded for j — 1 < /c < / + 1 and for / — 1 < /c < j + 1 
respectively by 

^-my\-e,-mw-e.^^ _ „ ^ fmax(2(-+^)^ 1) max(2(-+^)', 1), 

2 ^,sup5fc X |^^^(2("™+^)M)max(2(™+^)M), ^^"'^^ 

and since the exponential factors are bounded, after summation over / we get 

WR^gWrOL-- ^^sup^pigk, no.8) 



as desired. For R* in (10.1), we have 



< ^ ^m\j\+0j-m\k\-ek-m\l\-ei^j^_ 



k-l<j<l+l 
l-l<j<k+l 



;io.9) 



Then the sums over k and I are bounded in both cases by 

sup ^J, min(2™^-^^ 1) min(2-'"^■-^^ 1), (^g.lO) 

and hence \\R*f\\r~e-i];^^ < supif.(pigk, as desired. 
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Next we show the optimahty. Let 6 G Z, and choose any g which is piecewise 
constant on each dyadic interval (2-', 2-''^^), suppgf C p**, oo), and g >Q. Then for 
< r < 2^ we have 

/•oo na 

R^g{r) = h{r) I I h{s)~^(p{a)h{a)ag{s)dsda 
J2'' J2'' 

j-i (10.11) 

j>b k=b j>b 

where we denote g^ = ||5'||re-iL°°(r~2'=) ^ind (fj = ||v5||r~9Li(r~2J)- Choosing a test 
function i/.' G C^(0, 00) satisfying ip > 0, suppip C (0,2^) and {h \ ip) = 1, we 
see that (pj G i'^ {j > b) is necessary since we can choose arbitrary non-negative 
gk G £f (/c > b). Similarly by choosing suppf? C (0, 2''] and supp?/; C (2^^, 00), we see 
that (pj G (j < b) is also necessary. □ 

10.2. Double endpoint Strichartz estimate. Lemma 5.1 holds for more general 
radial potentials. We call 



-t 



^ ll^/IU?.. (10.12) 



t,X 



the Kato estimate for the operator H, and 

hh^L^) ^ I|m(0)||l2 + \\iut + HuWm^LVy (10.13) 

the double endpoint Strichartz estimate for H. Lemma 5.1 is a consequence of the 
following. 

Theorem 10.1. For any m > 0, the double endpoint Strichartz (10.13) holds for 
radially symmetric u{t,x) = u{t, \x\) and H = Ag"^"* = d^. + r~^dr — m^r'"^. 

Lemma 10.2. Suppose Hq and H = Hq + V are both self-adjoint on L^(M^) and 
\x\'^V{x) G L°°(M^). Assume that the Kato estimate (10.12) holds for H, and that the 
double endpoint Strichartz estimate (10.13) holds for Hq. Then we have the double 
endpoint Strichartz also for H. The same is true when we restrict all functions to 
radially symmetric ones, if V is also symmetric. 

Corollary 10.3. Let V = V{\x\) G C^(M^\{0}) be a radially- symmetric function 
with \x\'^V G L°°(M^), and suppose H = —A -\- V is self-adjoint on L^(M^). Let 
f{t,x) = fit, \x\) be radial. Then 

(1) the Kato estimate (10.12) holds for H if and only if the double- endpoint 
estimate (10.13) holds for IL , 

(2) both estimates hold provided 

infrV(r)>0, inf -r2(r1/(r))^ > 0. (10.14) 

r>0 r>0 

Since our linearized operator if* satisfies (10.14), the above implies Lemma 5.1. 

Proof of Corollary 10.3. The first statement follows directly from Theorem 10.1 and 
Lemma 10.2. For the second statement: the methods of [3], adapted to the 2- 
dimensional radial setting (detailed in [11]), imply that conditions (10.14) yield 
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the resolvent estimate (10.15), hence the Kato estimate, and the double endpoint 
estimate. □ 

Remark 1. While the double-endpoint estimate (10.13) always implies the Kato 
estimate (10.12), the reverse implication does not hold in general. For example, 
consider H = —A acting on 2D functions with zero angular average. The Kato 
estimate in this case can be verified, for example, by using the methods of [3] to 
establish the resolvent estimate 

sup||(i7-A)-VllL^.-<||¥'|l4,i, (10.15) 

from which the Kato estimate follows by Plancherel in t (see [11] for details). On the 
other hand, if the double-endpoint estimate were to hold for zero-angular-average 
functions, so would the endpoint homogeneous estimate. Since the latter is known 
to hold for radial functions (see Tao [18]), it would therefore hold for all 2D func- 
tions, which is false (see Montgomery-Smith [15], also see [18]). Alternatively, a 
constructive counter-example is given by placing delta functions of the same mass 
but opposite sign at (1, 0) and (0, 1) in the plane. 

Proof of Theorem 10.1. Following [14], we use the identity 

// F{s,t)dsdt = C dtF{s,t), 

JJs<t Jo ^ Jr ^ Ja-3r J a+r 

for the decomposition, where C > is some explicit positive constant. Define the 
bilinear operators /-,• for j G Z by 

a—r pa+Sr 

ds / dt(e^(*"^)^^'"V(s)k(t)>. 

— 3r J a+r 

where / and g are radial (i.e. /(s) = /(|a;|, s), etc.). 
The desired estimate follows from 

Using the 1/t decay for ||e 2 ||^i^j;^oo, we can easily bound the supremum of the 
summand. To get summability, we need decay both faster and slower than \jt. In 
fact we have, for >^ = </?(|x|) radial, 

lle^^^^^VlU-- < |tr^+1|^||^i,-M, -m < /i < 1/2. (10.16) 
This follows easily from the explicit fundamental solution 

(e-^"^V)(r) = f r e^^^'-^'y-J.^ (^) Mpdp 

J 

[cm a constant) in terms of the Bessel function of the first kind, for which 

sup s^Jm{s) < 00, — m < /i < 1/2. 

Next, when the decay is slower than 1/t, namely if we choose yU > in (10.16), 
then we get a non-endpoint Strichartz estimate using the Hardy-Littlewood-Sobolev 
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inequality in time: 

'"e-^'"V(.)rf.|U.< 11/11^,,^,..., (10.17) 

for < a < 1/2 and 1/p = 1/2 - a/2. 

Now the rest of the proof follows along the lines of Keel-Tao [13]. We will prove 
that 

\W^9)\ < 2^'^"^"^/1/llL?Li-<^ll^?llL|Li^-/', (10.18) 

for 

-m<a = p<0 (10.19) 

and for 

1/2. (10.20) 

For the first exponents (10.19), we use the decay estimate (10.16) and the — 
duality at each (s,t). Then we get 

J23 ' JR ' Ja-Zr J a+r F ~ I 



-I ~ [ ~'^^°'\\f\\L^{a-3r,a-r-,Ll--°')\\9\\L^{a+r,a+3r-,Ll'-'') 

23 ^ jR ^ 



- y2^°-|l/llL2,(_3,.<t_a<-r;L^'-")ll5'llL2,(r<t-a<3r;Li'-'^) 



where we used Holder for s, t, a. 

For the second exponents (10.20), we use the non-endpoint Strichartz (10.17) for 
both integrals in s and t, after applying the Schwartz inequality in x. Then we get 

\IjU.9)\ < / — / — ll/ILf(a_3r,a-r;Li--'')ll^?llLf(a+r,a+3r;Li'-) 



2J 



23 



r ./iB r 



L?(a-3r,a-r;L;^'"'')ll5'llL2(a+^^a+3r.il.-"), 



and the rest is the same as above, where = 1/2 + a/2 and 1/q' = 1/2 + /3/2. 
Thus we get (10.18) both for (10.19) and (10.20). By bilinear complex interpolation 
(cf. [2]), we can extend the region (a,/5) to the convex hull: 

1\ ^ / 1\ O 1 / N 

The only property we need is that this set includes a neighborhood of (0, 0), where 
we are looking for the summability. 

Now we use bilinear interpolation (see [2, Exercise 3.13.5(b)] and [16]) 

T:XiXX,^Y,+, (z,j,z+j G {0,1}) 

=^ T : Xe^^n X Xg^^r-i ^ ^ei+ea.n, V'^o = l/^^i + 1/^2, 
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where Xq^^ '■= (^o^^Oe.r denotes the real interpolation space. 
The above bound (10.18) can be written as 

\\I{f^9)\\i^i'-+P)/^ ^ \\f\\LlLl'~A\9\\L^^Ll--^^ 

where 1^ denotes the weighted space over Z: 

'■= l|2^"aj||^P(2)- 
Hence the bilinear interpolation implies that 

< II/IIl^l^.-II^IIl^lI'-"' (10.22) 

for all (a,/?) in (10.21), where 

and we used the interpolation property of weighted spaces (cf. [2]): 

(LP'", LP^P)e,g = a^p. 

By choosing a = (3 = in (10.22), we get the desired result. □ 

Proof of Lemma 10.2. By time translation, we can replace the interval of integration 
in (10.12) and (10.13) by (— oo,t). Then by taking the dual, we can also replace 
it by (t, oo). Adding those two, we can replace it by M. Then the standard TT* 
argument implies that 

l|e^''VllL?(L^,-i) < IblU^, l|e^^°VllL?(L^) < MW- 
Now let ^ 

u= f e'(*-^)^/(s)cis. 



Then the Duhamel formula for the equation 

iut + Hqu = f — Vu 

implies that 

^ = / e'^'-')^'^{f -Vu){s)ds. 



oo 



Applying (10.12) for H and (10.13) for Hq, and using L^'^ C L\, we get 

lkllL?(L-) < 11/ - yu\\Ll{L\) 

^ ll/IL?(L^) + lk^'^IU?=lkllL2(L2.-i) < ||/||l2(L2.1). 

Then by duality we also get 

Feeding this back into (10.23), we get 



(10.23) 



The estimate on e*^V is simpler, or can be derived from the above by the TT* 
argument. □ 
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Appendix A. Landau-Lifshitz maps from 

The same stability problem on instead of M^, is much easier in the dissipative 
case, because the eigenfunctions get additional decay from the curved metric on S^. 
Indeed we have convergence for all m > 1: 

Theorem A.l. Let m > 2, a G C and Rea > 0. Then there exists 6 > such that 
for any u{0,x) G with S{u{0)) < Am-K + 6^, we have a unique global solution 
■u G C([0, oo); Sm) satisfying Vu G L^;q^([0, oo); L"^). Moreover, for some fioo G C 
we have 

Mt) - e"^^^/I[/ioo]|Us° + £iu{t) - e'^'^fx^]) ^0 (t ^ oo). (A.l) 

Our proof does not give a uniform bound on 5 if m = 1, but we have 

Theorem A. 2. Let m = 1, a G C, Rea > and fiQ G C. Then there exists 6 > 
such that for any u{0,x) G Si with S{u{0) — /i[/io]) ^ we have a unique global 
solution u G C([0, oo); Si) satisfying Vu G ^^^^^([O, oo); L!^). Moreover, for some 
^oo E C we have 

\\u{t) - e'"^^/I[/ioo]|Ugo + g{il{t) - e'^'^^hM ^0 {t^ oo). (A.2) 

The proof is essentially a small subset of that in the case, so we just indicate 
necessary modifications. 

Outline of Proof. By the stereographic projection, we can translate the problem to 
with the metric g{x)dx^, where g{x) = (1 + r^/4)~^. The harmonic maps are the 
same, while the evolution equation is changed to 

qt = iSq — aL^g^^L*jq, S = g~^{q H 1^) ° iaL'tqdr. (A.3) 

Joo r 

In this setting we can use the "standard" orthogonality to decide 

= iz\ ghl), (A.4) 

since ghi G {r)^^L^. The energy identity 

dMh = -^a^ig^'Ltq\Llq) (A.5) 
implies the a priori bound on q: 

hhrL'i + \\g-'/'Ltq\\L2L^^ < WqmiLl ~ 6. (A.6) 
Since g~^^^ > r, we get (by using q = W*L^*q as on M^), 

MlI < WrLtqW^i < \\g~'/'Llq\\Li G L?. (A.7) 
Then by the orthogonality (A.4) we have z = R^^^L'^z with 

(t>s:=s'g{rs)h,/{ghl\hl), (A.8) 

and so 

\\z/r\\L2 < \\L'z\\l2Us\\lI- (A.9) 



32 STEPHEN GUSTAFSON, KENJI NAKANISHI, AND TAI-PENG TSAI 

Since 



min(l, s ^) (m > 2) 
min(s~"'^, (m = 1) 



g{rs)mm{r,l/r)-"'rdr ^ { . 1 li _2^ ) (A.IO) 



we have 



mk.~\' J"-'' (A.n) 

^ ' max s ,1 m = 1 



Anyway, if 6 is small enough (depending on s), we get by the same argument as on 



d2 



\\z\\x<\\q\\Lie^nLr. (A.12) 
Differentiating the orthogonality, we get 

KK I 9K) = -{MaLlq \ h\) - {gz \ {^rdr + tfi2K)h{) 

m 

= -{{Mr + '^M)aq\h\) (A.13) 
r 

- {gz I {^{rdr + m) + ifi^^hl - 
m 

where on the second equality we used that Uh\ = and {gz \ hi) = 0. Using that 

\{rdr + m)h\ + \{h- < min(r— i,r-3™-i) < (r)-^ (A.14) 

we can bound the last term in (A.13) by 

|/x|||z||l«> min(s^ 1), (A.15) 

which is much smaller than the term on the left. The second last term in (A.13) is 
bounded at each t by 

lklliil|/^?IU.^< Ikllii- (A.16) 

If m > 1, we can improve this for s < 1 as follows. By the same argument as on M^, 
we have 

Ml^^ < WmiLl < \\g-'^'Ltq\\L2, (A.17) 

where we need m > 1 for the boundedness of R^^ : r'^L^ — rL^ and i?^* : rLl L^. 
Then we can replace the above estimate in the region r < 1 by 

\\Q\\LfnLl\\K\\Ll^+L^ ^ Iklli2min(s^l). (A.18) 

Thus we obtain 

IIAIIlj < [""-^I, (A.19) 

* I C{s)d'^ (m = 1) 
and hence if 5 > is small enough, we get the desired convergence as on M^. □ 
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